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The experimental results for ZrZn2, URhGe, and in some pressure ranges also for UGe2, have 
shown that the ferromagnetic superconductors are weak itinerant ferromagnets. Guided by these 
results we describe the phenomenon of coexistence between equal spin triplet pairing superconduc- 
tivity (SC) and ferromagnetism (F) using the extended Stoner model, which includes in Hamiltonian 
^ , the on-site Coulomb interaction, U, and occupation dependent hopping integral. We use Hartree- 

, O ■ Fock (H-F) approximation and the Green functions technique. In the H-F approximation the on-site 

' Coulomb interaction plays the role of the on-site exchange (Hund's) field. All inter-site interactions 

will have included the inter-site kinetic correlation, {cf^Cja), within the H-F approximation. We 
introduce the pressure-dependence to the hopping integral. Numerical results are compared with 
experimental data for ZrZn2. The kinetic correlation creates the superconductivity without help 
of negative values of Coulomb interactions. The model can explain stimulation of triplet SC by 
the weak itinerant ferromagnetism. This effect was observed experimentally in ZrZn2. Numeri- 
5-H ' cal analysis also confirms the experimental effect of decrease in critical temperatures (Curie and 

lyj I superconducting) with increasing external pressure. 



PACS numbers: 71.10.Fd, 74.20.-z, 75.10.Lp 



Ch ' 1. INTRODUCTION 

o ■ 



The theoretical possibility of ferromagnetism (F) coexisting with the triplet parallel spins superconductivity (SC) 
was suggested by Fay and Appel for ZrZn2 [ij . Further theoretical development took place after finding experimental 
evidence for coexistence of triplet SC with F (see e.g. [1]). It is relatively recently that the so-called ferromagnetic 
QQ ' superconductors have been discovered, which at some highpressures exhibit the ferromagnetic and a spin triplet 
t — i; superconducting phases at the same time. At present UGe2 URhGe [3| and ZrZn2 [1] belong to the ferromagnetic 
superconductors. In ZrZn2 compound the ambient pressure strength affects the Curie temperature. This compound 
has the quasi-linear dependence of both magnetic moment and Curie temperature on pressure. 



00 . 2. THE MODEL HAMILTONIAN 

o 

^ ' Taking into account that the ferromagnetic superconductors are weak itinerant ferromagnets we can describe them 

|v>( ' by the extended Stoner model (see e.g. [6|), which includes the on-site Coulomb repulsion, U, and occupation dependent 

; I ' hopping integral. We can write the model Hamiltonian as 

^ \ ^ 4- 



<i3><y 



+ ^ ^nio-ni-o-, (1) 



where /i is the chemical potential. The At and t^x terms correspond to hopping and exchange-hopping interaction 
respectively, and are given by 

At - to - ti = to (1 - Si) , t,x = ^-^^ - ti = ito (1 + 51^2 - 2Si) , (2) 
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where Si = ti/to, and S2 = i2/ii- In equations above, to, ti, and t2 are the hopping ampUtudes for an electron with 
spin a when both sites i and j are empty, when one of the sites i or j is occupied by an electron with opposite spin, 
and when both sites i and j are occupied by electrons with opposite spin, respectively. 

In the calculations below, the constants to, Si and S2 will be assumed pressure-dependent. The kinetic interactions 
At and tex will also depend on pressure through Eq. @. The on-site Coulomb repulsion, U, will be assumed 
pressure- independent. Taking into account the results of [7| giving us the dependence of effective mass, to*, on the 
pressure, p, and comparing the dispersion relation in the tight binding approximation (at small k) with the expression 
£fe — h^k"^ l2m*, one obtains the relation 

^W-T^^ (3) 

where for ZrZn2 A = 0.017 ± 0.004 kbar"\ and the lattice constant a = 7.393A (see 3)- 

Further on the dependence 
of hopping integrals on pressure will be suppressed in the notation, i.e. t{p) = t, At{p) = At, tex{p) = tex- 

In Hamiltonian ([T]) there are terms with four and six operators. The terms with four operators will be approximated 
by the average of two of them multiplied by the remaining two. The averages of the spin-flip type, (cf^Cj-cr) , will be 
ignored. The six-operator term standing at tf,x will be approximated by the product of two averages of two operators 
multiplied by the remaining two operators (see Q). Using these approximations and dropping in the Hamiltonian 
spin singlet term and opposite spin triplet term we obtain, after transforming into momentum space, the following 
form 

H = Y,{el-^i + M-) n^, ^ (A^c+ cl,, + h.c.) , (4) 

k(7 k(7 

where = Suh^ is the spin-dependent modified dispersion relation, with the bandwidth factor, 6^, given by 

6" = 1 - Y [^tn^o - tex - - 2/,/_,)] . (5) 

The spin-dependent modified molecular field, Af^, is expressed as 

= Un_„ + 2z/_^ {At - 2texn„) , I„ = (c+ c^,,), (6) 

where z is the number of nearest-neighbors. The equal spin pairing (ESP) parameter, A^, for the two-dimensional 
square lattice is given by 

A^ = Atexl-a (A^ sin kx + a;; sin ky) , AJ(j^) = {c^+x(yyc^a) . (7) 

Solving the Green's function equations of motion (see [Q]) with Hamiltonian Q we obtain the following relation for 
the ESP superconducting critical temperature 

l=UexI-a—y\ ^tanh — ^ , (8) 



k 



where — \J (e^ — + Af'^)^ + (2A^)^, T is the temperature, and ks is the Boltzmann's constant. 
The equations for carrier concentration, n, and magnetization, to 



(9) 



can be obtained from n'^ 



n = — 
2 



N ^ 

k 



^ + ^:tanh.^^^ 



El 



2kBT 



(10) 



The Fock's parameter, I„ = 



{(^ta^j'y)' is given as 



2N ^ \zt 

k 



^^^^-^ + ^^%anh.^^^ 



E^k 



2fcsT 



(11) 
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3. NUMERICAL RESULTS 

Solving self-consistently Eqs ((8|)- pT|) we obtain the phase diagrams showing the dependence of superconducting 
critical temperature, ferromagnetic critical temperature, and magnetic moment on pressure. In numerical calculations 
we have used A = 0.013 kbar~^ in Eq. ([3]), and the following relations for 5*1 and iS'2 

_ ti ^ 0.35 _h ^ 0.25 

^ to l- 0.003p' ^ ti l- 0.003p' ^ ' 

The relation p2p was assumed as the relation similar to the relation describing the pressure dependence of the hop- 
ping integral t{p), Eq. ([3]). The value 0.003 in front of pressure was chosen as to obtain the Curie and superconducting 
critical temperatures in a good agreement with experimental data for ZrZn2. 

Fig.[l]shows the dependence of ESP superconducting, Tg^^ and ferromagnetic, T^" , critical temperatures on external 
pressure. The value of electron occupation n = 1.015 was chosen. This value allows obtaining the Curie temperature 
comparable with experimental data at p = kbar. The numerical results show that the Curie temperature at zero 
pressure is decreasing with growing carrier concentration (all remaining parameters being fixed). Note that the ESP 
superconducting critical temperature is magnified 10 times. As we can see, increasing the pressure causes decreasing 
of both superconducting and Curie temperatures. The pressure-dependence of critical temperatures is quasi- linear 
up to p « 16 kbar. The effect of two values of Curie temperature at a given pressure above 16 kbar is caused by the 
specific shape of the density of states (logarithmic) and by the H-F approximation used in the model. The proper 
curvature of the Curie temperature dependence on pressure can be obtained by assuming semi-elliptic density of states 
or by use of the higher-order approximations for four- and six-operator terms in the Hamiltonian (e.g. Hubbard III) 
(see (i). 

Fig. [2] shows the dependence of magnetic moment on external pressure. As we can see, the dependence rn(p) is also 
quasi-linear up to p « 16 kbar, in good agreement with experimental data. The strange behavior of this curve above 
16 kbar was explained earlier, in the description of Fig. [TJ 

CONCLUSIONS 

The simple itinerant electron model of Eq. JT]) with additional assumption given by Eq. ^ and Eq. (|12p can explain 
the pressure dependence of superconducting critical temperature, ferromagnetic critical temperature, and magnetic 
moment. 

The abrupt decrease of magnetic moment at p w 16.5 kbar is related to the structural phase transition (see [9!]), 
which we do not take into account in our model. 

It is worthwhile to note in here, that taking into account the inter-site correlations (given by the expression 
{cf^Cja)) gave us the nonzero values of parameters Si and 6*2, and equivalently the nonzero values of parameters At 
and tex- The last two are responsible for an existence of small superconductivity stimulated by ferromagnetism. Both 
superconductivity and ferromagnetism have pressure dependence in agreement with the experimental data. 
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FIG. 1: ESP superconducting (solid line) and ferromagnetic (dashed line) critical temperature vs. pressure. The following 
values were used: n — 1.015, U — 0.075 eV, and to — 0.125 eV. 




FIG. 2: Magnetic moment (in Bohr's magnetons) vs. pressure. The values used are the same as in Fig. (TJ i.e.: n = 1.015, 
U = 0.075 eV, and to = 0.125 eV. 



